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Abstract 

, In the phase space of the integrable Hamiltonian system with three degrees of freedom used to describe 

the motion of a Kowalevski- type top in a double constant force field, we point out the four-dimensional 
invariant manifold. It is shown that this manifold consists of critical motions generating a smooth sheet of 
, the bifurcation diagram, and the induced dynamic system is Hamiltonian with certain subset of points of 

degeneration of the symplectic structure. We find the transformation separating variables in this system. 
As a result, the solutions can be represented in terms of elliptic functions of time. The corresponding 
■ phase topology is completely described. 
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<^ Introduction 

, The equations of rotation of a rigid body about a fixed point in a double constant force field have the form 
00 ■ 

! I—— = la; X a; + ri X q; + r2 X /3, 

cn ■ (1) 

O ■ da dfB ^ 

> 



where ri and r2 are vectors immovable with respect to the body, a and /3 are vectors immovable in the 
inertial space, I is the tensor of inertia at the fixed point O, and u: is the instantaneous angular velocity (all 
' of these objects are expressed via their components relative to certain axes strictly attached to the body). 

It is assumed that the vectors ri and r-z have the origin O. The points specified by these vectors in 
the moving space are called the centers of rigging. 

System (1) is a Hamiltonian system in the phase space specified in R^{uj,a., f3) by geometric 
integrals; is diffeomorphic to the tangent bundle TSO{3). 

In [31, it is proposed to use the problem of motion of a magnetized rigid body in gravitational and 
magnetic fields and the problem of motion of a rigid body with constant distribution of electric charge in 
gravitational and electric fields as physical models of Eqs. (1). The results obtained for system (1) in [3] are 
presented in more details in 4 within the framework of investigation of the Euler equations on Lie algebras. 

In the general case ri x r2 7^ and a x /3 ^ 0, system (1) without additional restrictions imposed on 
the parameters, unlike the classical Euler-Poisson equations, is not reducible to a Hamiltonian system with 
two degrees of freedom and does not have known first integrals on P^ , except for the energy integral 

H = -Iu)-uj ~ ri-a r2'/3. 

In [Sj, the following assumptions are introduced for system (1): in the principal axes of the inertia 
tensor 

Oeie2e3, (2) 

the moments of inertia satisfy the conditions h = I2 = 2/3 and the vectors ri and Y2 are parallel to 
the equatorial plane Oeie2 and mutually orthogonal. For /3 = 0, the problem reduces to the Kowalevski 
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integrable case of rotation of a heavy rigid body [8]. Therefore, for the sake of brevity, the problem proposed 
in [3] is caUed the generahzed Kowalevski case and the problem with /3 = is called the classical Kowalevski 
case. 

By the appropriate choice of measurement units and axes (2) one can obtain 

I = diag{2,2,l}; (3) 

ri = ei, r2 = e2. (4) 

In [3], a new general integral is indicated for the generalized Kowalevski case. In virtue of relations (3) 
and (4), this integral admits a representation: 

K = {Lol-Lol+ai- [32f + {2uJiLJ2 +a2+ pif , (5) 

where uji,ai, and [3i {i — 1,2,3) are the components of the vectors u:,a, and (3 relative to the reference 
frame (2). 

In [l^, integral (5) is generalized to the case of motion of a gyrostat in a linear force field by sup- 
plementing a body having property (3) with an inner rotor generating a constant moment along the axis 
of dynamic symmetry. As shown, e.g, in [5], the component of the moment generated by potential forces 
introduced in [10] can be reduced to the same form as in Eqs. (1) with property (4). 

The complete Liouville integrability of the Kowalevski gyrostat in a double force field was proved 
in [2]. The Lax representation of equations of type (1) (with gyroscopic term in the moment of external 
forces, as in [10]) containing the spectral parameter was obtained under conditions (3) and (4). The spectral 
curve of this representation made it possible to find a new first integral which is in involution with the 
corresponding generalization of integral (5) and turns into the square of the momentum integral for (3 — 0. 
Multi-dimensional analogs of the Kowalevski problem were introduced in [2]. It was proposed to solve 
these problems by the method of finite-band integration. This program was realized in [2] for the classical 
Kowalevski top and new expressions for the phase variables in the form of special hyper-elliptic functions 
of time were obtained. The explicit integration of the problem of motion of the Kowalevski top in a double 
field and its qualitative and topological analysis have not been performed yet (see also a survey in [S]). 

The topological analysis of an integrable Hamiltonian system includes the description (in certain terms) 
of the foliation of its phase space into Liouville tori. In particular, this requires finding all separating cases. 
These cases correspond to the points of the bifurcation diagram of the integral map and, in the phase 
space, are formed by the trajectories completely consisting of the points at which the first integrals are not 
independent. 

In a system with three degrees of freedom, two-dimensional Liouville tori QifG, clS &- rule, filled with special 
motions corresponding to a point of the smooth two-dimensional sheet of the bifurcation diagram. The union 
of these tori over all points of the sheet is an invariant subset of the phase space. In the neighborhood of a 
point of general position, this subset is a four-dimensional manifold, and the dynamic system induced on this 
subset must be Hamiltonian with two degrees of freedom (degenerations of various kinds are expected at the 
boundary of this sheet or at the points of intersection of sheets) . Thus, the invariant subsets of maximum 
dimension formed by the points of dependence of integrals are specified (at least, locally) by systems of two 
invariant relations of the form 

/i = 0, /2 = 0. (6) 

The knowledge of all these systems and the analysis of the dynamics on invariant manifolds specified 
by these systems is essential to fulfil the topological analysis of the entire problem. 

In the generalized Kowalevski case, we know two systems of the form (6). The first one is obtained in 
[3]. Consider the manifold 

{K^OjcP^. (7) 

Due to the structure of function (5), this manifold is specified by two independent equations Zi = and 
Z2 = 0. An additional partial integral (Poisson bracket {^1,^2}) is indicated. The topological analysis of 
the induced Hamiltonian system with two degrees of freedom is carried out in [llj . It turns out that the 
invariant set is a four-dimensional manifold, which is smooth everywhere but the restriction to it of the 
symplectic structure degenerates on the set of zeros of the additional integral. This case generalizes the first 
Appelrot class (Delone class) [T of motions of the classical Kowalevski case. 
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The second system of the form (6) is obtained in [7j- It is shown that, for /3 = 0, the corresponding 
motions transform into so-caUed especially remarkable motions of the second and third Appehot classes. The 
present work is devoted to the investigation of the dynamic system on the invariant subset indicated in [7j. 

First, we make a general remark important for the sequel. The moment of external forces ri x Q: + r2 x /3 
appearing in (1) is preserved by the change 

where 6 is an arbitrary non-singular 2x2 matrix. Therefore, the a-priori assumption made in [5], [5] 
concerning the orthogonality of the radius vectors of the centers of rigging is redundant (it suffices to require 
that these centers lie in the equatorial plane of the body). This statement is trivial enough; the possibility 
of orthogonalization of any pair (ri,r2) or (a,/3) is indicated, e.g., in [5]. However, the authors of also 
indicate that, in general case, the second pair is not orthogonal. Moreover, in [31 HI HI [HI [7] , the angle 
between a and f3 remains arbitrary. This fact makes the corresponding formulas more complicated. 

Note that if the pair (ri, is made orthonormal, then there remains the arbitrary choice of G S0{2). 
Under such transformation, a new pair of radius vectors of the centers of rigging remains orthonormal and 
can be used as equatorial unit vectors of the principal axes (2) to preserve properties (4). At the same time, 
by choosing 

cos 6 sin 9 



e 



sin 6 cos 9 



tan 29= - 



we get the orthogonal pair (q:,/3). 

Thus, without loss of generality, in addition to relations (4), we can assume that the force fields are 
orthogonal. This simple fact has not been indicated yet. The elimination of the redundant parameter makes 
it possible to significantly simplify all subsequent calculations and to obtain results in a symmetric form. 

1 Invariant Subset and Its Properties 

In the sequel, we consider system (1) under assumptions (3) and (4) with the phase space specified by 
the formulas 

= a^, 13^ = 6^ a-f3 = 0. (9) 

The case a = 6 is singular. Indeed, as indicated in [lOj . in this case, there exist a group of symmetries 
generated by the transformations of the configuration space and, hence, a cyclic integral linear in the angular 
velocities. Therefore, for the sake of being definite, we set 

a>b. (10) 

Denote by G the general integral of the problem obtained in [2 and represent it in the form 

G^^{gl+gj) + ^uj39^-b^ai-a^f32, (11) 

where ga, gp, and g^y are the scalar products of the kinetic momentum lo; and the vectors ol, (3, a x /3, 
respectively. 

Introduce the function F by setting 

F = {2G-p^Hf -r'^K, 

where the parameters p and r are specified as follows: 

p^ = a^ + b\ r''=a^~ b\ (12) 

The latter is well defined due to inequality (10). Obviously, F is a combined first integral of Eqs. (1). 
Note that the zero level of the function F is specified by one of the conditions 

2G-p^H -r'^^/K = (13) 
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2G-p-'H + rWK = 0. (14) 

If /3 = 0, then these conditions reduce to the equations of the second and third Appekot classes, 
respectively [1]. 

Define the subset N C a.s the set of critical points of the function F lying on the level F = 0. 

The set N is definitely non-empty; it contains, e.g., all points of the form ui = UJ2 = 0, ai — P2 — 0, 
and a2 + /3i = 0, which are critical for K and turn the expression 2G — p^H into zero. 

The set N is invariant under the phase flow of (1) as the set of critical points of the general integral. 

The condition dF = means that the differentials of the functions H, K, and G are linearly dependent 
at the points of N. It immediately implies that the relation 

{2g-p^hf -r'^k^Q (15) 

for the constants of these integrals is the equation of one of the sheets of the bifurcation diagram (the 
investigation of this diagram has not been completed yet) of the generalized Kowalevski top. 

We use the following complex change of variables [7] (a generalization of the Kowalevski change; see 

my- _ 

= (ai - /32) + «(a2 + X2 = xi, 
yi = {ai+ /32) + i{a2~ (3i), y2^yi, ^^^^ 

Wl — LUl + iuj2, W2 = WT, W3 — 

Denote the operation of differentiation with respect to the imaginary time it by primes and rewrite the 
equations of motion in terms of the new variables: 

x[ = ~XiW3 + ZiWi, X2 ^ X2W3 - Z2W2, 

y'l = -yiW3 + Z2W1, y2 = y2W3 ~ Z1W2, ^^^^ 

2z{ ^ X1W2 - y2Wi, 2z'2 ^ -X2W1 + yiW2, 
2w[ = -{W1W3 + zi), 2w'2 = W2W3 + Z2, 2^3 = 2/2-2/1- 

Constraints (9) now take the form 

2 I 2 2 I 2 

21+2^1^2=?^, Z2+X2yi=r, 

9 (1°) 
X1X2 + yi2/2 + 2ziZ2 = 2p-. 

Further on, instead of integral (11), it is convenient to consider another general integral linearly ex- 
pressed via G and H, namely, 

M=^{2G-p^H). (19) 

On the level F = 0, we have 

K = r'^M^. (20) 
In terms of variables (16), rewrite the integrals H, K, and M as follows: 

H = wiW2 + ^wl-^{yi+y2), K ^ U1U2, 
M-~^F^ + ^iUi + U2), 
where ^ 

Fi = ^XiX2W3 ; ix2ZlWi + X1Z2W2), 

^JX\X2 

U^ = ^iwl+X,), U2^^{wi+X2). 
Xi X2 

Consider the function 

F2 = Ui- U2. 
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Proposition 1.1. In the domain X1X2 7^ 0, the invariant set N is specified by the following system of 
functionally independent equations: 

Fi =0, F2 = 0. (21) 

Proof. Represent relation (20) in the form 

[F^ - r\^,- Vlhf] [Ff - r^^. + ^^f] = 0, (22) 



where yUi and v C/2 are chosen to be complex conjugates of each other. 

On the level F = 0, the functions Fi, \/lJi, and \/c/2 are independent everywhere except for the set 

W1W2 = 0, xi = X2. (23) 

Therefore, the condition that the left-hand side of Eq. (22) has a critical point leads to Eqs. (21). It is 
clear that points (23) also satisfy these equations. Thus, it remains to notice that the functions Fi and F2 
are independent on the level F = Q everywhere in their domain of definition including points (23). □ 

The system of invariant relations (21) is obtained in |7j without using the first integrals. In virtue of 
the above definition and Proposition ll.il the indicated system describes a certain smooth four-dimensional 
(non-closed) manifold 

= NC^{xlX2 ^ 0}, 
and N is the least invariant subset of containing N^. 

Remark 1.1. It is easy to see that the invariant set iV, as a whole, is stratified, namely, 

N ^ t} N\ dim A^* = i, dN' C N' . 
i=i j=i 

Moreover, in virtue of Proposition [Ol all A^' with i < 4 belong to a subset of the phase space specified by 
the equation 

X1X2 = (24) 

(e.g., N'^ = {wiW2 = 0, KJ3 = 0,xiX2 = 0} is diffeomorphic to 25'"'^). Therefore, the existence of singularity 
(24) in the expressions for Fi and F2 is in no case accidental. If, in relations (21), we remove the denomina- 
tors, then the set of solutions of the obtained system contains the entire four-dimensional manifold (24). This 
manifold is not everywhere critical for the function F (however, F is identically zero on it). In particular, 
all trajectories starting from this manifold fill a set in P^, which is almost everywhere five-dimensional. 

The following statement demonstrates that if we restrict ourselves to relations (21), i.e., study the 
dynamics only on iV*, then we do not lose any trajectory of the dynamic system on N . 

Proposition 1.2. Set (24) does not contain subsets invariant under the phase flow of system (1). 

To prove this proposition, it is necessary to compute the derivatives of X1X2 in virtue of Eqs. (17) up 
to the fourth order, inclusively, and show that they cannot vanish simultaneously on the set specified by 
relation (24). It is worth noting that the indicated strong degeneration of this subset also takes place for 
motions of the heavy Kowalevski top. In that case, condition (24) means that the axis of dynamic symmetry 
of the top is vertical. Special attention is given to this phenomenon both in the classical papers (see, e.g., 
[1]) and in recent investigations dealing with the computer animation of motion (see [9] , where one can also 
find an extensive bibliography of works in this field devoted to the investigation of heavy Kowalevski tops). 

Proposition 1.3. The differential 2-form induced on the manifold by the symplectic structure of the 
space P^ providing the Hamiltonian property of Eqs. (1) is non- degenerate everywhere except for the subset 
specified by the equation L = 0, where 



1 

L = 



y/XiX2 



XiXr> + Z1Z2 

W1W2 H [Ui + U2) 
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Proof. The Poisson brackets of the functions on R^(u;,a,/3) specifying the indicated symplectic structure 
are determined according to the following rules [31 : 

{ai,aj} = {A,/3j} = {m,(3j} = 0, 

where gi = 2uji, g2 — 2u!2, and = 073 are the components of the kinetic momentum. 

In relations (25), we now pass to variables (16) and compute the Poisson bracket for the functions Fi 
and F2. In view of relations (21), this gives 

{Fi,F2} = -r''L. 

The tangent space TqN'^ is a skew-orthogonal complement of the span of vectors included in the 
Hamiltonian fields with Hamiltonians Fi and F2. By the Cartan formula (see, e.g., [6j, p. 231), the 
restriction of the symplectic structure to TqN'^ is non-degenerate provided that {Fi,F2}{q) ^ 0. 

□ 

Proposition 1.4. The function L is the first integral of the dynamic system induced on N"^. Moreover, this 
integral is in involution with the integral M . 

Proof. As shown in [7j, in virtue of (17) we can write 

Fi^fiiF2, F^ = fi2Fi. 

Here /ii and ^2 are functions smooth in the neighborhood of N*. In view of these equalities, apply the 
Jacobi identity to the functions H, Fi, and F2 and obtain that the double Poisson bracket {H, {Fi, F2}} is 
a linear combination of the functions Fi and F2. Therefore, L' = on the set specified by relation (21). 
It is shown by direct calculation that the following relation is true under conditions (21): 

= 2p^M^ + 2HM + 1 (26) 

and, therefore, 

L{L, M} = M{H, M} = 0. 

It means that {L, M} = for i 7^ 0. Hence, by continuity, {L, M} = everywhere on N'^. □ 

Thus, in the smooth part of the invariant subset N completely specifying the entire dynamics on iV, 
the equations of motion of the generalized Kowalevski top define the Hamiltonian system with two degrees 
of freedom with a closed subset of points of degeneration of the symplectic structure nowhere dense in N'^. 



2 Analytic Solution 

By Proposition II. 4i to integrate the equations of motion in the set N, we can use the integrals M and L. 
The original general integrals H, K, and G are expressed via these integrals by using relations (19), (20), 
and (26). 

Theorem 2.1. On an arbitrary integral manifold 

Jm,i = {M = m,L = e}(l N, (27) 

the equations of motion are separated in the variables 

X1X2 + Z1Z2 + r'^ X1X2 + Z1Z2 - r"^ 
si = ^ , , S2 = „ , (28) 

2^X1X2 2y/XiX2 

and take the form 



- ■ 's? - a2Wms2 ~ tsi + -^(P - 1), 

_< (29) 

V 4m 
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Proof. In virtue of the first equation in (21), the function M takes the following form on A'': 

In virtue of the second equation in (21), wc get Ui = U2- Therefore, the integral equation M = m 
yields 

Ui = r^m and U2 = r^m. (30) 
Determine Ws from the first equation in (21) and Wi and W2 from Eqs. (30). We obtain 

ziwi Z2W2 [x^ [¥[ 

W3= 1 , W2 = \ — Rl, Wl = \ — R2, (31) 

Xi X2 y Xi V ^2 

where 

Rl = \Jr'^m — x\ and R2 = \Jr'^m — X2- (32) 
Substituting these quantities in the equation of the integral L, we obtain 

rn(x^ + z^) — ix + V r^m? — 2r'^mx cos u + x'^ = 0. (33) 

The variables x^z, and a arc defined as follows 

x"^ ~ xiX2^ = Z1Z2, a:;i + X2 = 2a; cos (7, (34) 

and the radical in (33) corresponds to wxW2, and therefore is non-negative. The other radicals used above, 
including Ri and R2, are algebraic. 
Equation (33) now yields 

RiR2 =ix- m{x'^ + z"^), 



Rl + Rl = ^{[^a; - m{x^ + z"^)]^ - x^} + r^m, 



Introducing the polynomial 



$(s) = 4ms2 - Us +—(e'^- 1), 
m 



we can write in terms of variables (28) 



Ri + R2 = y/<^{s2) and Ri - R2 = v^$(si). (35) 

Si — 32 Si — S2 



Using constraints (18) and relations (34), we obtain 

1 

r 

Hence, in terms of variables (28), 



{zi ± Z2f = ^[{x^ + z^± r^f - 2x\p^ ± r% 



2r rz 2r 



zi + Z2 = Jsl-a^, zi-Z2 = ysi"^^- (36) 

Sl — S2 * Si — S2 * 

We now differentiate relations (28) in virtue of system (17). In view of (31), we get 

s'l = ^{zi + Z2){Ri - R2) and s'2 ^ ^{zi - Z2){Ri + R2) ■ 

Substituting expressions (35) and (36) in these equalities, we arrive at system (29). □ 

Remark 2.1. It is clear that the deduced equations can be integrated in elliptic fimctions of time. By 
using the standard procedure, the solutions are expressed in terms of Jacobi functions. Their specific form 
depends on the location of the roots of the polynomials under the radicals on the right-hand sides. The 
bifurcation solutions of systems of this type correspond to stationary points of one of the equations, i.e., to 
the cases for which the polynomial 

(s^ - a2)(s2 - h'^)^{s) (37) 
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possesses a multiple root. 

For dimension reasons, the original phase variables on manifold (27) are expressed via si and S2, though, 
in general, these expressions might be multi-valued functions. We now show that the latter have a fairly 
simple algebraic form. 

Introduce the following notation: 

Si = ^Jsl- a2, ipi = y-$(si), 

(38) 

5*2 = ^Jh^- si, = \fHs2); 

V' = 4msiS2-2^(si + S2) + ^(^^-l)- (39) 

Theorem 2.2. On the common level of the first integrals (27), hy using notation (38), (39), the phase 
variables of the generalized Kowalevski case can be expressed, in terms of variables (28), as follows: 

1 



2(Sl - S2) 

1 



012 = TTT 77[(S1S2 - a'^)(pi(p2 - S1S21P], 

^(Sl — S2) 

^1 = "■^T? ~ ^[(■SlS2 - b'^)(pi(p2 - SlS2'ip], 

Z[Si — S2) 

02 = ^ ,^ [{S1S2 - 6')V + SiS2ipiip2], (40) 

2(si - S2) 

r r 

0:3 = 5*1, 03 = S2, 

Si - S2 Si - S2 

r r 
wi = — -(£- 2msi)(^2, ^2 = 777 -(^- 2ms2)¥'i, 

2(si-S2) 2(si-S2) 

W3 = ^- (5'2¥'i - S^^p2)■ 

Sl - S2 

Proof. By using notation (12), we represent the compatibility conditions of constraints (18) in the variables 
X and z as follows: 

a;^ + + > 2a\x\ , \x^ + z^ - ^ 2b\x\ , 
whence we get the natoraZ ranges for variables (28): 

.s^ ^ and s^^b^. (41) 

Hence, rewriting Eqs. (29) in the real form, we conclude that, for given m and i, the domain of possible 
motions in the plane (si,S2) is determined, along with inequalities (41), by the inequalities 

$(si) < and $(52) > 0. (42) 

Thus, in particular, all values (38) are real on the trajectories of the analyzed system. The expressions 
for the complex variables xi, X2, z\, Z2, w\, W2, and in terms of si and S2 are obtained by application, 
in sequence, of relations (35) with regard for (32), then (36), and, finally, (31). After this, the variables yi 
and 2/2 are determined from the first two relations in (18). By the change of variables inverse to (16), we 
arrive at the required dependencies (40). □ 

Note that the value si = 00 is ordinary for the first equation in (29) (because the degree of the 

polynomial under the radical on the right-hand side is even). Moreover, if this value belongs to the domain 
of possible motions, then it is reached during a finite period of time. Similarly, relations (40) also do not 
have singularities in this case. This can be proved by the change of variables si 1— > 1/si. Thus, in particular, 
we have deduced analytic expressions for all cases in which the trajectories pass the set specified by relation 
(24). It means that we have constructed the complete analytic solution of the problem on the invariant set 
N. 
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3 Phase Topology 



In the regular case, the integral manifold Jm/ consists of two-dimensional Liouville tori. The cases when 
they topologically rearrange generate the bifurcation diagram of the system on N. It is natural to study this 
diagram in the plane of constants of the used integrals, i.e., as the set of critical values of the map 

J = Af X L : AT ^ R2. (43) 

Theorem 3.1. The bifurcation diagram of map (43) is a part of the system of straight lines 

e = -2ma ± 1, ^ = 2ma ±1, £^ -2mb ± 1, ^ = 2mb ± 1, (44) 

and the coordinate axes of the plane {£, m) lying in the half-plane i ^ and specified by the conditions of 
existence of real solutions 

£ ^ max{2ma — 1, —2mb + 1), to > 0; 

£s^-2to6+1, m<0; (45) 

£^1, TO = 0. 

Proof. According to Remark 12.11 the diagram contains the discriminant set of polynomial (37) formed by 
straight lines (44) (in the part corresponding to the existence of motions). 

The points of the coordinate axes in the plane {m,£) belonging to J{N) must be included in the 
diagram; indeed, it can be shown that the values m — and £ = are attained, in particular, on the subsets 
N'^ (i < 4), where N fails to be smooth (see Rcmark ll.ip . The analytical foundation for this inclusion is as 
follows. 

In Eqs. (29), we can pass to the limit as to ^ 0. As a result, by using relations (26), we conclude 
that 1^1 — > 1 and (£^ — l)/2m h. At the same time, the degree of the polynomials under the radicals 
decreases to three; the form of solutions changes. Moreover, it is clear that K equals zero on the set 
Nr\{M = 0}, i.e., the corresponding motions also belong to the class (7). [It is worth noting that, as shown 
in [llj . the restriction of symplectic structure to manifold (7) degenerates just at the points of corresponding 
trajectories.] Therefore, the value to = should be regarded as corresponding to bifurcation. 

On the other hand, the integral surface Jm.i is preserved by the inversion 

{as, (33,0:3) 1-^ (-aa, -/Js, -ujs). 

In relations (40), this inversion is realized either by changing the sign of the radicals Si and S2 or by 
the substitution (£,si,S2) '-^ {—£, —si, —32)- This means that Jm,i and Jm,-e are the same subset of the 
phase space. Therefore, we need to restrict ourselves to the values of £ of the same sign (to be definite, 
we choose non- negative values); the axis £ = becomes the outer boundary of the domain of existence of 
motions in the plane of constants of the integrals. In virtue of Eq. (33), £ can equal zero only for negative 
values of to. 

Thus, in addition to (44), the diagram should also be supplemented with the point {to = 0, £ = 1} and 
the semi-axis 

{£ = 0,TO<0}. (46) 

By analyzing the compatibility of conditions (41) and (42), we determine the actual region of existence 
of motions in the form (45). □ 

In Fig. 1, the domains with numbers 1-9 defined by the diagram in the plane {m,£) correspond to 
different types of the integral surfaces (27). The motion is impossible in the shaded region. 

To determine the number of tori for the regular manifold, we note that relations (40) give a one- valued 
dependence of the phase variables on two collections of quantities 

{si,Si,(pi) and {32,82, ip2)- (47) 

In this case, the signs of the radicals in (38) on each Jm,,i are arbitrary. However, along the trajectory, 
some radicals turn to zero and then change the sign. This means that two points that only differ by the sign 
of such radical lie in the same connected component of Jm,e- Therefore, the number of connected components 
of the regular integral manifold is equal to 2", where n is the number of quantities (38) non-zero along the 
trajectory. The value of n is determined according to the location of roots of polynomial (37) and does not 
exceed 2. 
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Figure 1: Bifurcation diagram and the domains of existence of motions 

Proposition 3.1. Assume that the analyzed domains are numbered as indicated in Fig. 1. Then the integral 
manifolds Jm.i can he described as follows: a) T^, in domains 1 and 8, b) 2T'^, in domains 2, 7, and 9, and 
c) 4T'^, in domains 3-6. 

To determine the type of critical integral surfaces, we note that, in each three-dimensional space of 
collections (47), equalities (38) specify a pair of cylinders (elliptic or hyperbolic) with mutually orthogonal 
generatrices. For the points of the straight lines (44), a pair of cylinders corresponding to one of the variables 
si or S2 has a tangency point. Hence, the line of their intersection is the eight curve \J S^. Thus, on 
segments of the straight lines (44) bounded by the points of lines intersection and internal for domain (45), 
the integral surface consists of the components homeomorphic to the product x {S^ V S^). Crossing such 
segment, we observe one of bifurcations 2T^ typical for systems with two degrees of freedom. The 

number of connected components of the form x (5^ V S*^) in the critical Jm,e is determined by the number 
of tori in the adjacent domains. Actually, the critical periodic trajectories (the traces of centers of the eight 
curve) are motions in which one of variables si or S2 remains constant and equal to the multiple root of the 
corresponding polynomial under the radical. In this case, either 5i = and = or 5*2 = and ip2 = 0. 
Hence, it follows from relations (40) that lu2 = co^ = in the first case and lui = uj^ = in the second 
case. The body performs pendulum motions in which the radius vector of one of the centers of rigging is 
permanently directed along the corresponding force field. In approaching the outer boundary of domain (45) 
with the exception of the half- line (46), the tori degenerate into circles (periodic solutions of the indicated 
pendulum type) and the surfaces x {S^ V S^) degenerate into eight curves. 

It is clear that the critical single- frequency motions do not appear in the half- line (46). The correspond- 
ing bifurcation in the segments adjacent to domains 5 and 6 is characterized by the fact that the number of 
connected components of Jm.o is half as large as at the close regular point of the plane (m, £). These are so- 
called minimal tori. The transition from domain 4 to a segment of the boundary set (46) is not accompanied 
by the decrease in the number of components of Jm.e f^nd all cycles homotopic to a certain marked cycle 
are folded so that each component covers the limiting component twice. In a sufficiently smooth case (e.g., 
in the case when L is a Bott integral on the corresponding smooth level of the integral M), a Klein bottle 
should be obtained as a result (see, e.g., |6J). However, according to the explicit equations (40), this is not 
true in our case. Most likely, this phenomenon is connected with the degeneration of the induced symplectic 
structure. 

Finally, consider the nodes denoted by Pi-Pa in Fig. 1. For these values of the constants of integrals, 
each surface Jm,t contains one singular point. These points correspond to the equilibria of the body in which 
both centers of rigging lie on the corresponding axes of force fields and, hence, the moment of forces is equal 
to zero. One of these points is stable: at Pi, the integral surface consists of a single point. The other three 
points are unstable. As indicated above, at the nodes P2 and P3 the integral surface is homeomorphic to an 
eight curve. At the node P4 the integral surface can be described as follows. Take a rectangle and identify its 
vertices with one point; it then can be filled with trajectories double-asymptotic to the singular point. The 
boundary of this set is a bunch of four circles. This boundary represents two pairs of pendulum motions; 
each pair is asymptotic to the highest position of one of the two centers of rigging. Take four copies of the 
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obtained set and attach the boundary of each to the same bunch of four circles. 

All indicated phenomena are readily established by analyzing relations (40) and the mutual location of 
the cylinders formed in spaces (47). 

Conclusions 

In the present work, we perform the complete investigation of motions of the generalized Kowalevski top 

playing the role of critical motions for the entire problem and generating bifurcations of three-dimensional 
Liouville tori along paths crossing the sheet specified by Eq. (15) of the bifurcation diagram S C R"^ of the 
general integrals of the problem. Inequalities (45) are used to deduce the equations of the boundary of a 
part of this sheet corresponding to the existence of actual critical motions, i.e., contained in E. 

Consider relation (22). It plays the role of the equation of the entire integral surface in the phase space 
for the collection of constants of integrals satisfying relation (15). It then follows, similar to the case of 
the heavy Kowalevski top (the second and third Appelrot classes), that the straight line {k = 0,2g = p^h} 
splits sheet (15) into two classes. In the first class specified by relation (13) and corresponding to the first 
non- negatively definite factor in (22), the obtained integral manifolds, being critical for the original system, 
exhaust the entire corresponding integral surface in as the limit of a concentric family of three-dimensional 
tori and are, in this sense, stable in P^. In the second class specified by relation (14) and corresponding 
to the second (hyperbolic) factor in (22), all obtained critical surfaces in P^ are hyperbolically unstable: 
on the same level of the first three integrals, one can find trajectories consisting of regular points and 
double-asymptotic to the corresponding two-dimensional tori of the system on the investigated invariant set. 
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